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Abstract 

We prove the existence of a 1/N expansion to all orders in /3 matrix models with 
a confining, offcritical potential corresponding to an equilibrium measure with a con- 
nected support. Thus, the coefficients of the expansion can be obtained recursively by 
the "topological recursion" derived in [CE06]. Our method relies on the combination 
of a priori bounds on the correlators and the study of Schwinger-Dyson equations, 
thanks to the uses of classical complex analysis techniques. These a priori bounds can 
be derived following [dMPS95, Joh98, KS10] or for strictly convex potentials by using 
concentration of measure [AGZ10, Section 2.3]. Doing so, we extend the strategy of 
[GMS07], from the hermitian models (/3 = 2) and perturbative potentials, to general 
(3 models. The existence of the first correction in 1/N was considered in [Joh98] and 
more recently in [KS10]. Here, by taking similar hypotheses, we extend the result to 
H all orders in 1/N. 
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1 Introduction 



1.1 Definitions 



We consider the probability measure /%« on Mr given by: 
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[&_, 6+] is an interval of the real line, — oo < 6- < 6+ < +oo, and /3 is a positive number. 
For (3 = 1,2,4, this is the measure induced on the eigenvalues of $ by the probability 
measure d§ e~~ TrV ^ on a vector space £^,/s of N x N matrices. Sn,i (resp. £n,2, 
and £n,a) is the space of real symmetric (resp. hermitian, and quaternionic self-dual) 
matrices [Meh04]. For general /3 > 0, when V is quadratic (Hermite weight), or log + 
linear (Laguerre weight), Dumitriu and Edelman have found [DE02] a measure on a set 
of tridiagonal matrices which induces the measure dfi N 'a on eigenvalues. This can be 
generalized to any even polynomial potential [Rid]. This was subsequently exploited to 
study these particular (3 matrix models by Ramirez, Rider and Virag [RRV06]. Though, 
there is no known plain random matrix whose spectrum is distributed according to 
d/j, N 'a for general p and V, we still speak of a "matrix model", and we call A; the 
eigenvalues. 

We define the unnormalized empirical measure M^ of the eigenvalues given by 
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and their Cauchy-Stieltjes transform, which are generating series for the moments of 
Mjv- In fact, we prefer to work with the generating series of the cumulants. They are 
also called "correlators" , and are defined for xi, . . . , x n 6 C\[fe_, b + ] by 



W^ b -M( Xl , 
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In particular, we have 
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V;[b-,b+] 
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When there is no confusion, we may omit to write the dependence in V and [&_, b + ] in 
the exponent. 

It is well-known, see e.g. [Dei99, Theorem6] or [AGZ10, Theorem 2.6.1 and Corol- 
lary 2.6.3], that 

Theorem 1.1 Assume that V : [&_,&+] — y R is a continuous function, and ifb T = too 
is infinite, assume that: 

hm inf — rVr > 1 

x^too 2 In \x\ 

If V depends on N , assume also that V — y V^ in the space of continuous function 
over [b—,b+] for the sup norm. Then, the normalized empirical measure Ljq = N~ l M^ 
converges almost surely and in expectation towards the unique probability measure 
/x e q := /^eq "' + on [&_,&+] which minimizes: 

e\p\ = J dpi(Ov {0 Ho - ff Mt)Mv) ln le - v\ 

Moreover, // eq has compact support. It is characterized by the existence of a constant 
C such that: 



Vx G [&_,&+] 2 /£ d/i cq (0 In |x - ^1 - ^ {0} (a;) < C , 

/i oq almost surely 2 J b + d/x e q(0 m |x — ^| — V^°j(x) = C . 

In particular, for any x G C \ [fe_, 6 + ], we have 

lim 1 W^(x) = / ^^ := ^ { - 1} (x) 
and the convergence is uniform in any compact of C \ [6_, 6 + ]. 

1.2 Main results 
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Our goal is to prove an asymptotic expansion in powers of 1/N when N — y oo for 

7 V;[b. 
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the partition function Z N 'l ' +J and the correlators W n ' (xi, ■ ■ ■ ,x n ). This is not 



always expected. In particular it is false when the support of /i eq "' , the limiting 
eigenvalue distribution, is not connected: corrections to the leading order feature a 
quasi periodic behavior with N (see [Eyn09] for a general heuristic argument). Our 
proof uses a priori bounds on the correlators, and what we really need is to establish 
that W n G 0(1) for n > 2. We shall prove this condition either based on a result of 
Boutet de Monvel, Pastur and Shcherbina [dMPS95] (also used recently in the context 
of ensembles by Kriecherbauer and Shcherbina [KS10]) or on a priori concentration 
of measure bounds when V is strictly convex. Our basic assumptions and main results 
are: 



Hypothesis 1.1 

• (Regularity) V : [6_,6 + ] — > M. is continuous, and if V depends on N, it has a 
limit V^ in the space of continuous functions over [b- , b + ] for the sup norm. 

• (Confinement) Ifb T = too, lim inf x ^ TOO ^^ > 1. 

• (One-cut regime) The support of /i eq ~' consists in a unique interval [a_, a + ] C 
[b-,b + ]- 

• (Control of large deviations) The function x G [&_, ft+]\]a_, a + [t— > \V(x) — 
Jin \x — ^|^ eq (0 achieves its minimum value at a_ and a + only. 



(Offcriticality) S(x) > whenever x G [a_,a+], where: 

S(x) = n 



d/X e q 



rir'sHard^-^r') 



dx 



UreSoft( x - Or) 



where r G Hard (resp. r G Soh) iff b T = a T (resp. r(b T — a T ) > 0). 

• (Analyticity) V can be extended as a holomorphic function in some open neigh- 
borhood of [a_, a+]. 

• V has a 1/N expansion in this neighborhood, in the sense of Hyp. 4-5. 

Notice that the "one-cut regime", "offcriticality" and "control of large deviations" 
assumptions automatically hold when V is strictly convex (see [Joh98, Proposition 
3.1], which extends easily to analytic functions instead of polynomials). 

Proposition 1.1 Assume Hyp. 1.1. Then, W n "' admits an asymptotic expansion 
when N — > oo: 

W^ b -M( Xl ,...,x n )= J2 N- k W^ k \ Xl ,...,x n ) 

k>n-2 

which has the precise meaning that, for all K > n — 2: 

A 



W^ b -M( Xl} ...,x n )= J2 N- k WV> {k} (xx,...,x n ) + o(N- 



k=n-2 

The o(N~ K ) is uniform for Xi, . . . ,x n in any compact o/(C\[6_, b + ]) n , but not uniform 
in n and K. Moreover, if (b T — a T )r > (meaning that a T is a soft edge), the functions 
W n are independent ofb T chosen such that (b T — a T )r > and Hypotheses 1.1 hold. 



Proposition 1.2 Assume Hyp. 1.1, and 6_ < a_ < a + < b + (all edges are soft). 
Then, Z N ^' + ' admits an asymptotic expansion when N — > oo; 



(^) ' exp E^'< ;W («) 



\fe>-2 



In oi/ier words: 

ViT > -2 Z v N f™ = Z N ^ E (^=^J exp £ N " K^ + o{N~ K 



k=-2 



Moreover, the coefficients Fa are independent of 6_ and b + chosen such that b- < 
a_ < a + < b + and Hypotheses 1.1 hold. 

Zn,g/3e is the partition function of the Gaussian (3 ensemble, defined by the 
quadratic potential Vq{x) = y. It is given by a Selberg integral [Sel44] (also in 
[Meh04]): 



" j >-•- ■< '-' • ^ rij=i r(i + j/3/2) 

•'— U '"' l "^ ; f (1 + p/2) N 



(2tt)^ 2 (^/2)-' w/ ^^ w 11j= V ' 7' ' (1-4) 



For hard edges (i.e. 6_ = a_ or b + = a+), one may still interpolate between Z N '^~' + ' 
and a Gaussian /3 ensemble restricted to some interval (Corollary 5.1), but the partition 
function of the latter is not a Selberg integral and thus not known in closed form. 

1.2.1 Commentary 

When V does not depend on N and p, W n has a very simple dependence in j3: 

[k/2\+l . v l_ s _ n 9 N fc+2-2a-n 

^: ;W = E (§) (!-|) wr^ fe+2 - 2 ^ (i-s) 

and likewise: 

Lfc/2j + l . !_ fl 9Nfc+ 2- 2o 

F v ;{fc} = ^, + x: (0 (i - 1) ^^2-*, (1 . 6) 

Assuming existence of the 1/iV expansion, or at the level of formal matrix integrals, 
the recursive computation of the W n and J^'fow was developed by Chekhov and 
Eynard in [CE06]. For /3 = 2, it is well-known that Eqn. 1-3 is an expansion in even 
powers of N, i.e. Fgj} 2 ' — 0. Such a result goes back to the so-called topological 
expansion of t'Hooft, shown in the context of matrix models by Brezin, Itzykson, 



Parisi and Zuber [BIPZ78]. Indeed when /3 = 2, the sum in Eqn. 1-5 has only one 
term, namely k = 1g — 2 + n, which is present only when k = n mod 2, and likewise 
for Eqn. 1-6 which can be considered as the case n = 0. 

At the asymptotic level, the case /3 = 2 was tackled in [APS01]. For /3 = 1,2,4, 
the partition function and the correlators can be computed with the help of or- 
thogonal polynomials [Meh04]. These are solutions of a Riemann-Hilbert prob- 
lem [FIK92], for which the large N asymptotics have been intensively studied 
[BI99, DKM+97, DKM+99b, DKM+99a] with the steepest descent method introduced 
in [DZ95]. As a consequence, Ercolani and McLaughlin [EM03] were able to prove the 
existence of a 1/N 2 expansion of In Zjj B=2 . However, the topological expansions in the 
cases /3 = 1 and 4 are technically more involved in this framework, and have resisted 
to analysis up to now. Let us mention however that universality of the fluctuations of 
the eigenvalues could be obtained in these cases, see [DG07, DG09]. 

Integrability properties of /3 matrix models are unraveled for general (3 > 0, in par- 
ticular there is no known orthogonal polynomials techniques to evaluate the partition 
function Zjjo and the correlation functions W n (xi, . . . ,x n ). Yet, it is always possible 
to study the Cauchy-Stieltjes transform of the empirical measure of the eigenvalues 
and the "loop equations", also called Schwinger-Dyson equations or Pastur equations 
[Pas 72], that govern its expectations and cumulants. Thanks to the rough bounds for 
WY and W 2 V established in [dMPS95], Johansson [Joh98] proved a central limit theo- 
rem and obtained the first correction to W\ when V is an even polynomial satisfying 
Hyp. 1.1. This was also the subject of a recent work by Kriecherbauer and Shcherbina 
[KS10], with Hyp. 1.1 only. These authors have obtained in particular the expansion 
of InZ^ B up to a 0(1) when N — > oo (see their Theorem 2). 

The determination of H/f ;{_1} [Wig58, BIPZ78, AG97] and W 2 ' {0} [AM90, Bee94, 
Joh98] has been known for long, in /3 ensembles or many other matrix models. It 
was also observed long ago [AM90] that, if a 1/N expansion is assumed to exist, the 
loop equations turn into a system of recursive linear equations determining fully the 
decaying orders. To solve it, one just has to invert a linear operator /C. Recursiveness 
is a consequence of the assumption or the fact that W n G 0(N 2 ~ n ), which allows the 
determination of the leading order of W n -\ without knowledge of W n (for n > 3). These 
techniques found their origin in [AM90, ACM92, ACKM93, ACKM95] and culminated 
with the formalism of the "topological recursion" of [Eyn04, EO07] for = 2, and 
[CE06] for any fixed (3 > 0. 

In this article, we observe that /C _1 is a continuous operator on some appropriate 
space of analytic functions. Combining with the a priori control on correlators which 
dates back to [dMPS95], we prove the existence of the full expansion. 



For strictly convex potentials, concentration inequalities also provide rough bounds 
on the correlators. In this framework, loop equations were used in [GMS07] to establish 
the asymptotic expansion of models of several hermitian random matrices (/3 = 2) 
with strictly convex interactions. Maurel-Segala [MS] also studied models of several 
symmetric random matrices (/3 = 1) with strictly convex interactions. In order to prove 
the asymptotic expansion, the main step of [GMS07] was to show that some operator 
on non-commutative polynomials could be inverted, with bounded appropriate norm, 
and this was only done in a perturbative regime. Here, thanks to complex analysis, 
the potential need not be a small perturbation of the quadratic potential. 

Our techniques could also be applied to other matrix models. For instance, the 
convergent /3, 0(n) matrix model: 



N 



d ^,AO(n)W- 7 V;R + ll dA * e n M , \ W2 

Z 7V,/3,0(n) i=l LLl<i,j<N\ A ^ A 3J 

An important point is that the corresponding quadratic functional: 

£\p] = If MOdp(v) \ - § in |e - v\ + £ i" ie + v\] + 1 / MO v {0} (0 

J Jr 2 + vi i j i j R+ 

is strictly convex in the regime |n| < /3, therefore ensuring unicity of its minimizer. 
Besides, the analytic tools for the recursive determination of the one-cut solution to 
the loop equations of the 0(n) model in this regime were clarified in [BE11]. The 
existence of a full 1/iV expansion for convergent 0(n) matrix models under Hyp. 1.1 
could probably be established by following the lines we are presenting for the (3 matrix 
models. 

An open challenge, which would be interesting for further applications, is to extend 
the same reasoning to complex measures, i.e. to Eqn. 1-1 where the potential V is 
complex-valued, or/and where the eigenvalues are integrated on a given path in the 
complex plane. 

Outline of the article 

We first study in Section 2 the weak dependence in the bounds of integration under 
weak assumptions on V. In particular, we may trade the initial interval [&_,&+] for a 
finite interval [o_,a+]. We then write in Section 3 the corresponding loop equations 
for the correlators. Section 4 is devoted to the proof of the asymptotic expansion of 
the correlators with slightly stronger assumptions (Prop. 4.1). They are weakened in 
Section 5 to complete the proof of our main results for the correlators (Prop. 1.1) and 
the free energy (Prop. 1.2). We also remind how early steps of our proof imply the 
central limit theorem of Johansson (Prop. 5.2). 
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2 Weak dependence on the soft edges 

In this section we show that the partition function and the correlators depend very 
weakly on the boundary points of the interval of integration [6_, b+] if they are soft, 
i.e. do not coincide with the boundary points of the support [a_,a+] of the equilib- 
rium measure. We show more precisely that this dependence yields only exponentially 
small corrections, by deriving a large deviation principle for the law of the extreme 
eigenvalues. This point was already studied in [AGZ10, section 2.6.2] under a technical 
assumption [AGZ10, Assumption 2.6.5] that we replace here by assuming that the rate 
function of our large deviation principle vanishes only at a_ and a + . 

2.1 Large deviation principle for the extreme eigenvalues 

Under the assumptions of Theorem 1.1 on an interval [a_, a+], we can define: 

jV;[a.,a +]{x) = ^M _ j** ^"-^(fl In \ X - £\ 

when x G [a_,o+], and +oo otherwise. Suppose that [o_,a+] ^ [a_,a+], and set: 

jV;[a_,a + ]^ = J V;{a„,a + ]^_ ^ jV;{a„,a + ]^ 

t£[a-,a+] 

We define also Jmax~ (x) (resp. i7 mi „ ~' a+ (x)) which is equal to J v ^ a - ,a+ \x), except 
when x g] — oo, at-] (resp. [a+, +oo[) where we set its value to +oo. 

Proposition 2.1 Let V : [a_,a.f] — > M. be a continuous function, and if a T = too, 
assume that: 

liminfZ^) >i (2-1) 

x->too 2 In |x| 

Assume that t J v >l a -' a +] vanishes only at a_ and a + . Then: 

(i) /3 l 7max~' a+ (resp. /3^7 m S~ ) ^ s a 9°°d rate function on [a^,a + ], which vanishes 
at a + (resp. a-). 

(ii) The law of A max (resp. X m i n ) under fi^p satisfies a large deviation principle 

with speed N and rate function equal to f3Jm£k ~ (resp. ^ mi „ ' j on [a_, a + ]. 
In other words, for any closed subset F, or open subset fl, of [a_, a + ]: 

limsup 1 ln/i^- a+] (A max G F) < -0 MJ^-^(x) 



AT-»oc 



N '' xgf 



liminflln/i^' a+] (A max Gfi) > -0 inf J^- a + 
and similar statements hold for A min . 



In particular, for any e > 0, 



limsup — ln^pl" ,a+] (A min < a_ - e) < (2-2) 

limsup — ln/xj^"'^ 1 (A max > a+ + e) < (2-3) 



Proof. 



jV;[a^,a + ] j s fl g 00( ^ ra £ e f unc tion. J' v '<\- a -> a +} is lower semicontinuous as a supre- 
mum of the continuous functions 

Jp [a -' a+] (x) ■■= ^r~ f + d ^ la -' a+] (0 ln [max(|x-e|,e)]- c inf ^-^(0 

Moreover, by the assumption of Eqn. 2-1, it goes to infinity at infinity. Hence, 
jV;[a-,a + ] ^gg compact level sets. Since it is non-negative, it is a good rate 
function. 

The law of the extreme eigenvalues is exponentially tight, that is: 

limsup limsup — ln/4/t"' a+] (A max > M or A min < -M) = -oo (2-4) 

By [AGZ10, Lemma 2.6.7], it is enough to show that: 

^V;[a_,a + ] 

lim sup — ln .. , — r < oo (2-5) 

For this purpose, observe that by Jensen's inequality 

yV;[a-,a + ] f „„ : jV-1 ^ AT „ JV-1 

*N,P ,,V;[a_,a + ] 

7 V;[a-,a+] ~ ^iV-1,/3 






f + dA w exp (/3 ]T ln \X N - A 4 | - ^^(Ajv) - £ £ F(A 4 ) 

AT-l AT-1 

2 ^ ln |Ajv - Ai| - (JV - l)V(X N ) - £ F(A*) 



j = l 2=1 

where we denoted x the law on A^ given by: 



dx(x 



= l [a _, a+] (x)dx e _ fm k= f a+ ^ e ^y i0 



K -'a_ 



The function £ H y J R dx(Ajv) ln | X N — £| is bounded on compact sets and going to 
infinity like ln |£|, so is bounded from below, by a constant ^-. We can rewrite: 



yV;[a-,a+] 



J N,/3 



V;[a-,a + ], 






By exponential tightness [AGZ10, Eqn. 2.6.21], we know that there exists a 
constant k 2 > so that 

- &s +] [L»-i<y)] > -&t ] [ l ^(\v\)]>-^ 

So, if we set n 3 = x\Y] an d choose k 2 large enough, we have: 

yV;[a-,a+] 

-^ -> Ke -^ N ^ 8 

7 V;\a-,a+\ — 

with a positive constant 5 = —K\ + k^ + k 3 . This justifies Eqn. 2-5 and completes 
the proof of Eqn. 2-4. 

Upper bound for large deviation of the extreme eigenvalues. We give the argument 
for the minimal eigenvalue, the case of the maximal eigenvalue being similar. By 
exponential tightness (Eqn. 2-4), it is enough to prove a weak large deviation 
upper bound, that is control the probability of small balls. First, observe that 
for any x — a_ > 2e > 0, 



/% 



N,fi ' + [^min > X] < fl N ' t ' + [^jv(l[a_,a_+e]) — 0] 



is of order e Kt for some K e > by the large deviation principle for the law of 
Ln under fi N a~, see e.g. [AG97] or [AGZ10, Theorem 2.6.1]. Moreover, the 
probability that A min is smaller than a_ vanishes and therefore we have 

limsuplimsup — \n[i N o~' a+ (A min g] — oo, a_ — e] U [a_ + e, +00 [) = —00 

Hence, we may and shall concentrate on probability of deviating on [a_,a_], 
and actually we may restrict ourselves to the case where a_ and a + are finite by 
Eqn. 2-4. We let F be a closed subset of [a_, aJ\. We then have: 

^ [ ;-' a+] [\^ EF}=Y N [ dee-S^O EiV (£) 

Jf 



where we introduced: 



Z 



T V;[a-,a + ] 



v J N-l,/3 

Yn — 



7 V;[a-,a+] 
N,0 



AT-1 

En® ^/-^/SM^I-f«)JJ 1[fl _ A](6 



j=i 



The first step is to bound the term ^n{0 from above. Notice that the logarithm 
is uniformly bounded from above on compacts so that the exponent is at most of 
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j£rV;[a_,a+] 



order N. Therefore, we may and shall assume that under fi N Zi g , L^-i is 

at a distance smaller than k > from the equilibrium measure /i eq := jU e q , 

since the opposite event has probability smaller than e - r "-( N ~ 1 ) for some T K > 0, 
see e.g. [AGZ10, Theorem 2.6.1]. Thus, we have for large N: 

E N b(£,) < e~ T * N2/2 + e ^ N -^ BU Pd^^)< K (-^F+Jinie-^ld^qW) 

where we take the supremum over probability measures on [a_,o + ]. We observe 
also that for all probability measures \x on [a_, a + ], and for any ( > 0: 

J ]n\Z-7i\dfjL(Ti)<<f> c (jjL,€)= ln[max(|^-7;|,C)]d/*(77) 

where <f>((fjL, f ) is continuous in \x and £, and <j>^(fj, eq , f ) converges towards 0o(/^eq, 
as ( goes to zero. We deduce that: 

limsupsup sup /3 (/ In |f - 7]\dnirj) - ^- j < -/3 inf j^; [«-."+](£) 
Therefore, for any 7/ > 0, and iV large enough, we conclude that: 



(r/~Pmi ieF J v ^- a +\z)) 



supE^f) < e ^'-^^ "' +J ^ (2 _ 6) 

The second step is to bound Yn from above. We observe that, for any e > 
small enough, and any x G [a_ + e,a + — e], there exists <5 £ going to zero with £ 
so that 

7 V;[n-,a + ] 
1 __ Z 7V,/3 






( .s i "- wi (r* c ^"«'ni£-^ 



,,- + : W-l 






'/%-l„3 



where we have finally used Jensen's inequality. But A — > ^ J ^ha |f — A|df is 
bounded continuous on [o_, a + ] and therefore by the large deviation principle for 
the law of the empirical measure Lat_i under fi^Z\ 3 (with rate function 

which vanishes only at /i eq ) we deduce that: 

> lee 



1 - • -efv(x)-2N6e p (tf-l)/&(jT£ln|e-A|dO<WA) 
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Hence, by taking e sufficiently small independently of N, and optimizing over the 
choice of x e]a_, a+[, we conclude that for any 77" > 0, and N large enough, 



1 . -TV 

> e 



(y'+/3inf es[a _ ia+] J V ^- a +\0) (2-7) 



Putting Eqn. 2-6 and 2-7 together, we deduce that for all 5 > and N large 
enough: 

V;[a-,a+] /^ g ^ < giV/3 (-inW 1 7 Vi[o -«H- ] (*)+to'ee[a_,« + ] J v * a -< a +Ht)+s) 

which provides the announced upper bound. As a consequence, since we assumed 
that the rate function only vanishes at a_, a + we deduce that for any e > 0, there 
exists 5 t > so that: 

^S~ ,a+] (Amin<a--e)<e-^ (2-8) 

as well as a similar result for the largest eigenvalue. 

Lower bound for large deviation of extreme eigenvalues To establish a lower 
bound, we start again from Eqn. 4-11 with an open ball B =}x—e, x+e[c [a_, a_]: 



V,[a-,a+] 



(Kin eB)=Y N [ dCe-i^S^O 
Jb 



but replace the role of Y N and E N in the bounds. Namely, we first have by 
Jensen's inequality: 

in0w w (tf > /fArP /dx«.A)(i8 5:£l 1 li'lC-A J |-f(JV-l)V(0) 



dee-5^^S2v(e)>«Jve 



with 



Kjv = / d£e ^^<_i,/3 [!A min >d 

Thanks to Eqn. 2-8 (note that it applies similarly to NV/(N — 1) as the assump- 
tions does not depend on the fine asymptotics of V), we know that n^ converges 
towards a non vanishing constant. Moreover, the logarithm, once integrated 
against d£, produces a smooth bounded function and therefore we can use the 
convergence of Lat_i towards /x eq under /i^Ij a to conclude that: 

hminf-ln / d£ e -f v «>5tf > ~ V . ~ , ym 
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Letting now e going to zero in B =]x — e, x + e[ proves that: 
liminfliminf — In /" d^e - ^ Sjv(0 > fi( [ d/i cq (ri) \n \£ - r]\ - —^-) 

e ^ N^oo N J B (x, e ) W 2 / 

(2-9) 
To bound Y/v from below, it is enough to bound 1/Y/v from above, which can be 
done in the same way we bounded S^r from above in the argument for the upper 
bound. We finally conclude: 

liminfliminf — ln/^-' a+] (A min e]x -e,x + e[)> -f3J v > [a -> a+] (x) 
which completes the proof of the large deviation principle. o 



2.2 Weak dependence on the soft edges 

We first state the global version of our result: 

Proposition 2.2 Let V : [a_,a + ] — > R be a continuous function, and if a T = too, 
assume that: 

hm inf - . . > 1 

x^roo 2 In \x\ 

Suppose a_ < a_ ; and assume furthermore that the minimum value of l J v 'l a -' a +] is 
achieved only on [a_,a + ]. Then, for any e > 0, there exists rj e > so that: 

4J- a+] = 4J- £ ' a+ l (1 + 0(e- N *)) , 

and there exists a universal constant 7„ > such that, for any x\,...,x n G 
(C\[a_,a + ])»: 

\W^ a -^( Xl , . . . ,x n ) - W^-^( Xl , ...,x n )\< ln \, n ( 2 " 10 ) 

[U =1 d{xi,[a_,a + \) 

A similar result holds for the upper edge. 
We also have a local version: 

Proposition 2.3 Let V : [a_,a + ] — > R be a continuous function, and if a T = too, 
assume that: 

hm inf - . . > 1 

x^tod 2 In \x\ 

Suppose a_ < a + , and assume furthermore that the minimum value of J v >\ a -' a +\ is 
achieved only on [a_, a+\. For any e > small enough, there exists n £ > so that, for 
any a'_ e]a_, a_ — e[: 



o i r 7 V;[a'_,a+] 

d a >_ In Z N 
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< e ~ Nrk 



and, for any Xi, . . . , x n G (C \ [a'_, a+\) 
Wx 1 ,...,x n eC\ [a'_,a+], 



d a ' W n '" " + (Xi,...,X n ) 



< 



JnN r 



-Nr, e 



Y\ n i=1 d(x il [a'_ 1 a + ]) 
A similar statement holds for derivatives with respect to the upper bound. 

Proof. If a_ 7^ a_, let a'_ G]a_,a_[. Notice that: 

7 V;[a'_,a+]' 



N,P 



,V;[a~,a + ] 



V;[a_,a + ]r x / 

JV 8 L^min J^ W— 



/%,/3 



(2-H) 



N,$ 



If now : [o_, a+] — >■ C is a bounded continuous function, we can write: 



V;[a-,a+], 



V;[a'_,a+) 



7 V;[a'_,a+] 

J N,8 



V;[a'_,a + ] l 



Thus, we find: 

kS-' a+] [0(A)]-^ a '-' a+1 [0(A)]| < 2(sup Ae[a _ ia+]J v|0(A)|)^J- a+] [A mm <a'J 
This can be applied for the disconnected correlators: 



n N 



w 



[Xi, . . . ,x n ) — fl N 'p [ I I / v \ 

„■ i „■ 1 .7 2 



and we obtain: 



\W V J a -' a+ \x 1 ,...,x n )-W V n ' [a '-' a+ \ 



X\, . . . , X n J 2; 



j=l j,=i 



2N r 



W = id(xj,[a-,a + ]) 



V;[a_,a+]r x ^ / 

N 8 l^min _ u — 



/%,/3 



Similarly, one finds: 

la TT / -V;[a'_ l a + ] 

\Oa> Wn 



(r r)\< 2Nn 8, \nZ V '' [a '-' a+] 



The correlators W n ' a+ are just sums of monomials of the form 
Wni a+ (h) ■ ■ ■ Wn m a+ \l m ) where h,...,I m is a partition of {xi, . . . ,x n }. 
So, it is enough to establish the weak dependence at the level of the partition function. 
The global version is a direct consequence of Eqn. 2-2 applied to Eqn. 2-11: 



, z V;[a'_,a+] 

lim sup — In 1 ^ =■ 

-AT \ z V;[a-,a+] 



<0 



JV-+00 



y V,/3 
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For the local version, we rather need to bound: 



d a ' \nZ 



V;[a'_,a+] Ar ^N-1,/3 



gZr;\a'_,a+] 



N,p 



N V;[a'_,a + ] P'n-W 



P - 



JVV -^+Ef = 7 1 m|A i -a'J 



&[<*!, «M 



If a'_ G]a_,a_[ is fixed, by the large deviation principle for L N _ 1 under //^l\' 
since the logarithm is a lower semi continuous function, there exists 7 > such that, 
for any e > 0, for N large enough: 



/%-l,/3 



-^^ + Ef = - 1 1 ln|a'_-A,| 



<7e" 



-/3AT(l-e)J rl/; l a -'' I +)(a / _) 



Moreover, we have seen in Eqn. 2-7 that for iV large enough: 






7 V;[a'_,a+] 

J N,/3 



< ^(HWfel..,^] J^'- a+ \0 



By assumption, J^f-'^Vj = J^ a -' a+] (a'_) - inf ee[ a_,c*_] ^-^(O > 0, leading 
to: 



lim sup — In 

JV->oo N 



d n! InZ 



V;[a'_,a + ] 
N,/3 



<0 



which is the bound we sought. The arguments at the upper edge are similar. o 

3 Loop equations 

We shall assume in this Section and also in Section 4: 

Hypothesis 3.1 —00 < a_ < a + < +00. 

Indeed, considering finite intervals [a_,a.j-] is convenient to ensure from the beginning 
that the Cauchy-Stieltjes transform yields functions which are holomorphic outside 
[o_,a+]. We also assume in this section: 

Hypothesis 3.2 V : [a_,a + ] — > C can be extended as a holomorphic function in some 
open neighborhood of [a-,a + ]. 

This will allow us to use complex analysis (Cauchy residue formula, moving the con- 
tours, etc.) 

We shall derive the "loop equations", also called Schwinger-Dyson equations or 
Pastur equations [Pas72] in this context. These equations express the invariance by 
change of variable of an integration, up to boundary terms. We stress that these 
equations are exact for finite N. Although the technique is well-known, we recall the 
derivation here for the matrix models with edges a-,a + in order to have a self- 
contained presentation. 
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3.1 First version 

Theorem 3.1 Loop equation at rank 1. For any iGC\ [o_, a + ]: 

W 2 (x,x) + {W^x)) 2 + (l - l)^( w i« 

(x-a-){x-a+) \Jc([a-,a+]) 2in x ~£ {x-a-){x-a+) 1 J 

C([o_,a + ]) is a contour surrounding [a_,a + ] in positive orientation, and included in 
the domain where V is holomorphic. 

Theorem 3.2 Loop equation at rank n. Let xj = (xi)i e i a (n — l)-uple of spectator 
variables in (C \ [a_, a+]) n . For any igC\ [a_, a+]: 

W n+1 (x,x,xi) + ^2W\j\ +1 {x,xj)W n „\j\{x,x I \ J ) + (1- -)-—{W n (x,xi 



JCI 



_ N n p_ i (,- a _)(,- a+) 

\/c([a_,o + ]) 2nr x - £ (x - a_)(x - o+) 

,2^ d / ^- 1 (x > x AW )- ^n:ijr^ ) ) w w - 1 (x J ) , = o 

3 ^— ' dxj V x — Xj / 

Proof of Theorem 3.1. For any smooth real- valued function h, and e > small 
enough, 

Vv : A ^ A + e/i(A) 

defines a differentiable family of diffeomorphisms from [a_,a + ] to some interval 
iph >e ([a—,a+]). We assume hereafter that h(a-) = h(a + ) = so that ?/^ e ([a_, a+]) = 
[a_,a + ] for £ small enough. We have: 

1=/ d/^fl(V>v(Ai ),..., VM^at)) (3-1) 

When e — > 0, the first subleading order of the right hand side must vanish. It can 
be computed in three parts. A first term comes from the variation of the Lebesgue 
measure Y\ i dAj, which is given by the Jacobian of the change of variable: 

N N 

n^VMAi)) = (n dA *)(l + ^ / h'(Z)dM N (Z) + o(e] 

A second term comes from the variation of the Vandermonde: 

|A(^ £ (A))|^ = |A(A)|^l + ./3 £ h{K l I ^ + o(e) 

l<i<j<N % ■? 

= \A(X)f {i + e g ( ^ fe ® ~ h v iV) dM N (®dM N ( V ) - J h'(0dM N (0) + o(e) 
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The last term comes from the variation of the Boltzmann weight: 

N N 



[]e^«)l=(n e -T^)(l-^ / V>(Uh(UdM N (i)-rot:: 



i=l i=l 

Summing all terms up, the first order in e in Eqn. 3-1 vanishes iff: 



V;[a-,a+] 



h{£) - h( v ) 



t-i) 



dM N (Z)dM N (rj) - N / V'(0h(Z)dM N (Z) 



1- 






V;[a_,a + ] 



ti(Z)dM N (£) 



(3-2) 



Note that even though this equation was obtained for real-valued functions h, we 
can at this point remove this condition by linearity. To obtain an equation involving 
correlators, one can take for x G C \ [o_ , a+] the function h defined by: 

h{£) = 7. = * h a_ + a + - x - 4 

x-£ x-f 

thus preserving [a_, a + ]. We recall that V is holomorphic in a neighborhood of [a_, a+]. 
So, by Cauchy formula, for any contour C([a_,o+]) surrounding [a_, o+] inside this 
neighborhood and not enclosing x: 



V'(Z)(Z-a-)(Z-a + ) 

x — £ 

Hence, we obtain: 



dMjv(0 = / dM N (£) 



W 2 (x,x) + (W 1 (x)) 

-N 



C([a-,a+}) 



N 2 



drj V'(t))(t) — Q>-){t) — a>+) 



(x — a_)(x — a+) 

d£ 1 (e-o(e-a+) 

(;i[ „ „ |D 2ivrx-£ (x-a_)(x-a+) 



^'(O^i (0 



H)(-sw*»- 



iV 



/37 V dx v ' y/ (x — a_)(x — a + ). 
Proof of Theorem 3.2. By definition of the cumulants, if we define a shifted potential 

2 



W n (x,x 2 ,...,x n ) 



I3N 



d e (w^r\x 2 ,... 



,i )■ 



e=0 



Notice that the matrix integral with this shifted potential is still convergent, because 
the eigenvalues live on the finite interval [a_,o+]. Therefore, we can obtain the loop 
equations at rank n by taking a perturbed potential in Thm. 3.1: 



V( Xr , t2) ,..., {Xn ; tn) (0=V(0+Yl 



and identifying the term in 



lli=2 [bn) 6 



i=2 

when €{ — ¥ 



Xi-i 
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3.2 Second version 

Here is another equivalent form of the loop equations. All W n depend implicitly on the 
interval of integration [a_, a + ]. 

Theorem 3.3 Loop equation at rank 1. For any x G C \ [a_, a+\: 

W 2 (x,x) + (W^x)) 2 + (l - ^j^(w x (x] 



-N 



c'([a_,a+]) 2«VT X - £ 



2 /3 a _lnZ <9 a+ lnZ 



/3 V x — a_ 



x — a i 







C([a_,a + ]) zs a contour surrounding [a_,a + ] in positive orientation, and included in 
the domain where V is holomorphic. 

Theorem 3.4 Loop equation at rank n. Let xj = (xj)j e / a (n — l)-uple of spectator 
variables in (C \ [a_, a + ]) n_1 . For any x G C \ [a_, a+]: 

W n+ i(x,x,x 7 ) +^Wjj|+i(a:,xj) W„_|j|(x,a;j\j) + ( 1 - -J -p \W n (x,xi) 

d£ V(0^„(^i 



JC7 



-AT 



C([a-,a+]) 



2Z7T 



X-£ 



+ 2 V— ( W ™-i( x ' x i\{i}) - w n-i{xi)\ 2 (d a _W n „ 1 (x I ) t d a+ W n -i{xi) 






ft 



+ 



x — a. 



x — a + 



Proof In the former proof, if we use a change of variable h which does not preserve 
[a_, a+], the partition function becomes (to first order in e): 



Z iV "^ "% 



J. + e ( fc(a_) 3 a _ \nZl' [a -' a+] + h(a+) d a+ InZ 



V;[a_,a+1\ ,- , 



Thus, Eqn. 3-2 receives those extra terms, and becomes: 



V;[a-,a+] 



HO - Kn) 



£-v 



K) 



2\ V;[a_,a+] 



dM N (OdM N (n) - N I V'(€)h(S)dM N (£) 

2 



ti(OdM N (0 +- Q [h{a„)d a ^Z V / a -' a+] ' *•'- ^ a i„ ^;[a-,a + . 



+ /i(a + )5 fc lnZ Ar 
In particular, when we choose /i(£) = ^37, we obtain: 

, 2 „ / de ho^mo 



W 2 (x,x) + (Wi(x)) -N 



l ~l)lw' 



e([a-,a+]) 2i7T Z - £ 

2 <9 a _ In Z 2 d a+ In Z 



/3 x — a_ /3 x — a + 
The loop equation at higher rank can be deduced as before by perturbing the potential. 



o 
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3.3 Remark 

If we compare those expressions to the first version of the loop equations, we find by 
consistency: 

and for higher correlators d aT W n -.i(xi) equals 



a- T -a T 2 J C ({a_,a + ])2w 



(f " a-. T )V'(£)W n (t) + E ^T ( (x * ~ a -rWn-i{xi)) I 



for r G {±}. 



4 The I/AT expansion 

4.1 Notations, assumptions, proposition 

This section relies on complex analysis and inequalities for probability measures. We 
make four assumptions on the potential V, which are valid only in this section. The 
link with our main theorem will be done in Section 5. 
We keep on with the assumption: 

Hypothesis 4.1 — oo < a_ < a + < +oo. 

Since V is smooth, the equilibrium measure djUeq ' (£) will in fact be a density 
p(x)dx, where p : [a_, a + ] — > [0, +oo] is a continuous function. We call supp L = {x E 
[a_,a + ] p(x) > 0} its support. In the hermitian case (/3 = 2), a 1/N expansion is 
expected only when supp L is connected. We assume here also: 

Hypothesis 4.2 V leads to a one-cut regime, i.e. the support of p e q is an in- 

terval [a_,a+] C [o_,a + ]. 

In order to write the loop equations as in Section 3, we assume: 

Hypothesis 4.3 V is real-valued on [a_,a + ], and can be extended as a holomorphic 
function on some open neighborhood U of [a_,a + ]. 

We justify in Remark 4.1 later that there exists a unique analytic function y : U — >■ 
C U {oo} such that, for any x G [a_, a+], we have: 

p(x) = — lim y(x + ie) (4-1) 

in e->o+ 



This function can be written y(x) = S(x)a(x), where S is now a holomorphic function 
defined on U, and a is of the form: 



^=1/ "^ aT \ (4-2) 



The lower edge a_ is 



• either a hard edge, meaning that a_ = a_. Then, p(x) G 0(fx — a_) ^ 2 ) when 
x —)• a_. 

• or a soft edge, meaning that a_ < a_. Then, p(x) G 0((x — a-) 1 / 2 ) when 
x — >■ a_. 

and the same distinction exists independently for the upper edge a + . Our discussion 
holds for both hard and soft cases. However, a key technical assumption is: 

Hypothesis 4.4 V is off critical on [a_,a + ], in the sense that S(x) remains positive 
on [o_, a+]. 

For instance, Hyp. 4.2 and 4.4 automatically hold when V^ is strictly convex. For a 
generic V satisfying Hyp. 4.2, we have S(a-) > and S(a + ) > 0, so we can always find 
an interval [a_,a-).] which is a strict enlargment of [a_,a+], such that Hyp. 4.4 holds 
on [a_, a + ]. We call "critical point on [o_, a+]", the situation corresponding to a choice 
of V such that S has a zero on [o_, a + ]. In this article, we do not tackle the question 
of the double scaling limit for /3 matrix models (N —¥ +00 and coefficients of V finely 
tuned with A^ to achieve a critical point when iV = 00). Though, this would be a very 
interesting regime in relation with universality questions, considering the absence of 
Riemann-Hilbert techniques when (3 7^ 1,2,4. 

We shall allow V itself to depend on N and have a 1/N expansion. To give precise 
statements about those expansions, we need some notations. For any Jordan curve T, 
we note Ext(r) (resp. Int(r)) the unbounded (resp. bounded) connected component 
of C \ T. In the following, we fix once for all a Jordan curve Te, and a sequence of 
nested Jordan curves (Ti)i eN , which all live in C \ [a_, a + ], and such that 

(i) T E C U. 

{11) {xeu S(x) = 0} n int(r E ) = 0. 
(Hi) vign r, c int(r m ). 
(iv) vign r ; c int(r B ). 
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The contour configuration is depicted in Fig. 1, where the zeroes of S were called S{. In 
the remaining of the text, F will refer to a Jordan curve in Int(r^)\[a_, a + ]. We will use 
the following norm on the space % n \ a a i of holomorphic functions on (C \ [a_, a+]) n , 
which behave as 0{\jxi) when X{ — > oo. 

Il/llr= SU P |/(ari,...,x n )| = sup|/(xi,...,x„)| 

a;jeExt(r) x t er 

The second equality is a consequence of the maximum principle. One can easily derive 
the following useful inequalities: 



tW 



V/ettJ^] V/GN 



/(•) - /(xo) 



-So 



< HA™ < G 



(4-3) 



where ^ = 27rd2( - r r — y is a finite constant depending only on the relative position of 
Ti and T m . 



c 


X 

• 


Si 

• 


/ 

/ 

f 
1 
J 

/ singularities of V f 

X 
X 




— ""^"T ' \ 








"^ r? ^i±lrC\\ 


V Q 




:= ======^))) ) 




~ 




" 


V 


• 


/ X 

U = domain of analyticity of V 



Figure 1: Hypothesis on the location of s; and contour configurations 



Now, we can state our last assumption: 

Hypothesis 4.5 V admits a 1/N asymptotic expansion: 



v(x) = Y, N ~ kv{k} ( 

1 E 

k>0 



X) 
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with functions I/W independent of N, such that yW is holomorphic in U . The =r E 
equality means that, for any positive integer K , there exists a positive constant vk such 
that, for N large enough: 

K 



sup 

«er B 



v K 



v(o - J2 N ~ k y {k Ho < n- {k+i) 

(The maximum principle implies automatically the same bound with Y replacing Ye)- 

In many applications, V is independent of N (i.e. V = V^). There is however no 
difficulty in our reasoning to consider potentials which depend on N within Hyp. 4.5. 

Our intermediate result is: 

Theorem 4.1 If Hyp. ^. 1-^.5 hold, the correlators admit an asymptotic expansion 
when N — > oo with respect to the norm || ■ || r , of the form: 



Vn > 1, W n = J2 N ~ k W n } 



k>n-2 



where W^ } en^ a _ >a+] . 

4.2 Relevant linear operators 

4.2.1 The operator /C 

(2) 

We introduce the following linear operator defined on the space ti^.L a i of holomor- 
phic functions on C \ [o_, a + ] which behave as 0(1 /x 2 ) when x —¥ oo: 

(Kf)(x) = 2W { 1 - 1 \x)f(x)-l ^M(-J_ + c W°>)'(£)/(fl 

Jc(la-,a+]) 2«vr L(x) \x-£ J 

This operator for an appropriate choice of L and c appears in the loop equations. We 
have found the following choice convenient: 

L(x) = I [ (x — a T ) 

reHard 

if Soft = {±} or Hard = {±} 

— - — if r G Soft and (-t) G Hard 



We may also rewrite: 



(£/)(*) = -2y(x) f{x) + ^j^ (4-4) 



with 

(e/)(ar) = - 

>C([a-,a + ])UC(x) 



(Q f) (x) = - I H L(0 (-J- + c) (V^)'(O /(£) 
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where C(x) is a contour surrounding x only (computing a residue at x). For any 
/ G KL a _ i, (Q f) is analytic, with singularities only where (V^)' has singularities, 
in particular is holomorphic in the neighborhood of [a_, a + ]. We have set: 



y {x) = -Wl~ Li (x) + 



{-i>,^ , (^ { ° } )'(^ 



2 

y is discontinuous on the support of /i eq (see Thm. 1.1), i.e. on [a_, a+] C [a_, a+], 

but analytic on C\[a_, aq_]. We justify in Remark 4.1 that y(x) = S(x)a(x) where a(x) 
was introduced in Eqn. 4-2 and the squareroot is chosen with its usual discontinuity 
on M_. Let us call s, ^ a_, a + the zeroes of S(x) in the complex plane, and we assume 
that they do not lie in [a_,o + ] (Hyp. 4.4). 

It is clear that Im/C C %\.L_ a y Here, W^ (hence y) has only cut [a_, a+], and 
this operator is invertible 1 Its inverse can be explicitly written, it is given by Tricomi 
formula [Tri57]: 

VxG<C\[a_,a + ], VsGlm/C, {KT x g){x)=l ^- ^— ^f\ |^r 

Jc([a.,a+]) 2t7T £-X <j{x) 2j/(f) 

(4-5) 

where <x(x) = y(x — ct_)(x — a+), and where we integrate over a contour surrounding 
[o_, o+] but not x. Indeed, if g G Im/C, we can write for any x G C \ [a_, a+]: 

?(*)/(*) = Res-^a(0/(0 

d£ ?(Q /(Q 

C([o_,a+]) 2l7T £ - X 

d£ g(Q 1 ( a(n , (g/)(0" 
C([a_,o f ])2i7re-a;23/(0 ^ y ^ } L(0 

d£ g(g 9(0 

C([a.,a+]) 2m i-X 2j/(£) 

In the second line, we moved the contour from a neighborhood of x to a neighborhood 
of [a_, a + ], and used the fact that a(£) G 0(£) and /(£) G 0(l/£ 2 ), so that the residue 
at oo of the integrand vanishes. In the fourth line, we use the fact that L is chosen 
such that (flue) — sbj, which is holomorphic in a neighborhood of [o_, a+] thanks to 
Hyp. 4.4. Since, (Qf) is also holomorphic in a neighborhood of [a_,a.f], the contour 
integral of this term vanishes. For our purposes, it is not necessary to describe the 
vector space Im/C. Notice that if we apply /C _1 to a function g G Im/C which is 
furthermore holomorphic outside C \ [a_, a+], we can contract the contour C([a_, a + }) 
to a contour C([a_, «+]). 



In general, on the space of holomorphic functions with g + 1 cuts, dimKer/C = g, and one has to 
prescribe g cycle integrals in order to define an inverse operator. 
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4.2.2 Continuity of K and K' 1 

The key fact in this article is that /C _1 is a continuous operator in (im/C, || • || r ): 

Lemma 4.1 Im/C is closed subspace ofH^! a , for the topology induced by the norm 
|| ■ || r , and there exists a constant k > 0, such that: 

YgGlm/C, WIC' 1 g\\ r < k\\g\\ r 

We call ||/C _1 || r; the infimum of such constants k. 

Proof. Let us prove first that /C, as a endomorphism of Ti^.L a i, is continuous. 

For any / G Ti^.L a i in formula 4-4, if x runs along T, we can move the contour 
C([a-, a+\) U C(x) to T^ and find the bound: 

V/G||/C/|| r < 2(max ge r| y (x)|)||/|| r + ^ ma ^ er ^^ 1 (^-pr + c) max|(yW/(6| 

27T mm :re r|L(x)| \d(r£,r) / £er B 



< 



2m^ 6r | 9W | + Si) -ggg woi ( i + x {0) , ■ 

27r nun lS r \L(x)\ \d{T,r E ) > £er B 



(4-6) 



We have used again the maximum principle for g to find the second line. Likewise, we 
can show that /C~ x : Im/C — > H; , , is continuous. In formula 4-5, we put x on 
T, and move the contour from C([a_,a + ]) to T E in Eqn. 4-5. Doing so, we pick up a 
simple pole at £ = x, and we find: 

, r -i v x iM. 1 / de 1 £(flg(fl 
^ ^ W ~ 2y(x) + 3F(x) 7 r£ 2*vr £ - x 25(0 

We find the bound: 

nr-i I, ^ llffllr Z( t e) max eerB |L(£)| ll5llr B 

11^" ffllr — r> .__;._ I../_M "T 



2 min^gr |y(x)| 47rd(r,r B ) min xer |cr(x)| min ?e r B |S(£)| 

\2 min^gr | y (a;) | 47rd(r,r s ) min^r |cr(x)| min ?er \S({,)\) r 

(4-7) 

where we used the maximum principle in the last line. Eventually, let us show that 
Im/C is a closed subspace of fi\_.L_ a i- We pick up a sequence (g n ) n in Im/C converging 
towards g G n{.L_ a i for a norm || • || r on a given contour r . Let (f n ) n be a sequence 
in % V L_ a i such that g n = )Cf n , or equivalently f n = IC~ 1 g n . Using Eqn. 4-7 for any 
contour T, we know that ||/ n || r < k ||<7n(|r f° r some constant k > 0. So, f n is a locally 
bounded subsequence of holomorphic functions in C \ [a_,o+]. By Montel's theorem, 
it admits a subsequence (/^,( n ))n converging to some / G 7^./ a i uniformly on any 
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compact of C\ [o_, a+]. Then using Eqn. 4-6, a^( n ) = JCf^n) — > /C/ for the norm || • || r . 
In particular, g(x) = K,f{x) for all x G Ext(ro). Since g and / are both analytic in 
C \ [a_, a + ], they must coincide on C \ [a_, a + ]. Hence, g G Im/C, showing that Im/C 
is closed. o 

||/C _1 || r is controlled by the distance of the zeroes Sj to the support [o_,a+]. This 
motivates Hyp. 4.4. 



4.2.3 The endomorphism "negative part" 

Let g be a holomorphic function at least in a neighborhood of [o_,a+]. The following 

) 

[o_,a+] 



endomorphism of %\.L_ a i often appears in the loop equations: 



We may write sometimes Af g [f(x)] as an abuse of notation. The analyticity assumption 
on g ensures that M g is a continuous operator with respect to the norm || • || r . Indeed, 
let us put x on Y and move the contour C([a^,a + }) to T E : 

KUM = »<»)/(*) + £ § f| (^ + =) 9 (0 /«) 

Thus, the maximum principle implies: 

Uf(f)\\ < IMIIIfll + i ^-( 1 + ld > l maX ^ rJL(e)l llqll 

WMgWWr < ll<7llrll/llr + 2yr U(r B ,r) + |C 7 min xer |L(x)| lM ^ 



< 



I II i ^V X i I A max ge r B |L(Q I .. .. 

^ Ilr + 2tt U(r*,r)' MC V min xer |L(a;)| ^ llr * 



4.3 Order of magnitude of W n 

If there exists a 1/JV expansion, W n ought to be of order of magnitude N 2 ~ n . Let 
us start with a lemma explaining how this can be infered from rough bounds on W n . 
Hereafter, Oi(- ■ ■ ) or o/(- • • ) mean 0(- ■ ■ ) or o(- ■ ■ ) with respect to the norm || • || r . 
Since the contours T; are ordered from the interior to the exterior, being a oi+\(- ■ ■ ) is 
weaker than being a oj(- • - ). When the index / is not precised, it is understood that 
the bound holds for any integer I. 

Lemma 4.2 Let <LiWi := N~ l Wx - w}~ 1} and I > 0. Assume <LiWi G oj(l), and 

/or a// integer n > 2, assume W n G Oi(N). Then: 



Vn>2 ||W n || r4n _ 6+i G 0(iV 



■2-n"\ 
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Proof. Let 5qV = V — V™. Let us rewrite the first version of the loop equation at 
rank n > 2 (Thm. 3.2): 






W n {x, x T ) = A n+1 + B n + C n _! + D n _ 



where: 



C*n-1 



-^ oVy [/(x)]+2(5_ 1 H/ 1 )(x)/(x) 



--^W„+i(a;,a;,x/) 



1 

'iV 



E 



rtl,n 2 >2 JC//|J|=ni 
"1 +n2=n+l 



12 E" 



N p ^ dxi 



^ W / ni (x,x J )W / n2 (x,x/\j) 
7|J|=m 

W n _i(x,a;7\ W ) 



1 2 

iV/3 



E 

-eSoft 



^Wn-lfal) 



L[Xi 






+ c)W n -i{x I ) 



X 



Firstly, as <l_iWi and (5oV)' goes to uniformly on T_i when N — > oo, we observe 
that for any fixed integer k, and A" large enough: 



[i-e N , k+1 )\\w n \\ Th+i < 



\KT 



- fc+i 



/C + 5/C + 



1 

iV 



2\ d 

1 ~^)dx 



W n 



- fc+i 



+ 



N 



G \\W, 



n\\r k 



(4-8) 



where 



£jV,Jfc+l 



\K~ 



\r k+1 (W(s vy\\ rk+1 +2\\5^W 1 \\ rk+1 ) 

goes to zero as iV goes to infinity for k + 1 > /by assumption. G is defined in 
Eqn. 4-3. We assume hereafter that iV is large enough so that EN,k+i is smaller than 
1/2. Secondly, we know from Proposition 2.3 that -D n -i £ 0(e~ Nv ), so this term does 
not contribute at any order of magnitude N~ k . Now, if we assume that W n G Oi(N) 
for all n > 2 (this is obviously true for n = 1), we always have A n+ i G 0/(1) and 
C n _i G 0/ +2 (l), whereas the last last term in Eqn. 4-8 is bounded by hypothesis for 
k>l. 

Now, we want to bound W n by induction on n. At rank n = 2, we have B 2 = 0, 
and we deduce from Eqn. 4-8 that W 2 G 0/ +2 (l). Then at rank n — 3, the product 
term 5 3 is 0, +2 (l/iV) and C 2 is 0, +4 (l/JV), thus W3 G /+4 (A 4 ) = O e (l). Then 
similarly at rank n = 4, the product term £> 4 is 1+4(1 /N) and C3 is Oi+q(1/N), thus 
W 4 G O j+8 (1). This implies in return that A 4 G Oi+&(l/N), thus W 3 G O i+6 (l/iV). 
And so on . . . The result can be proved by a triangular induction, as depicted in Fig. 2. 



25 



w 2 
w 3 

W A 
W 6 

w 7 
w 8 



0{N) 




\ 



0(1/N) 

OWN) 

i 
0(1/N) 

i 
0(1/N) 

i 
0(1/N) 



0(1/N 2 ) 

I 
0(1/N 2 ) 

0(1/N 2 ) 



Figure 2: Triangular recursion for Lemma 4.2 



At each vertical step, we are forced to trade the contour T^ with the exterior contour 
Tk+2 in order to control the C terms. So, to go from W n G Ok n {N 2 ~ n ) (in the n th 
column) to W n G Ok n+1 (N 2 ~( n+1,> ) (in the (n + l) th column), we must reach W n+ 2 in 
the n th column. This is done by two vertical steps, thus k n+ i = k n + 4. Since k<i = 1 + 2, 
we have k n = An — 6 + / for all n > 2. o 

Lemma 4.3 // there exists 7 G [0, 1[ and 5 G [0, oof such that W n G O (N"* n ~ 5 ) for all 
n > 2, then: 

\\W n \\ T e 0{N 2 ~ n ) 

wherel = 2\(-f- 1 - l)" 1 ] . 

Proof. Now, let us rather assume the existence of 7 G [0, 1[ and 5 > such that, 
for all n > 2, W n G Oo(N in ~ s ). D n _\ being always exponentially small, it does 
not matter in our discussion. At rank n = 2, as B2 = and C\ G 02(1), we have 
W 2 G 2 (max[-^ W 3 , 1]). We also have for all n: 

A n+l e O (m&x[N^- 5 -^\l}) 

B n G O (max[A^ 7n - 2<5 -( 1 ^ 7 ), 1]) 
C n _x G 2 (max[i\H™-Mi+7) ; x]) 
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When these 0(- ■ ■ ) decay, it does not hurt to consider them as 0(1). So, our bounds 
are upgraded at least to W n G 2 (max[iV 7n ~' 5 , 1]) with 8' = 5 + 1 — 7 > 5. By repeating 
the argument k times, we obtain for all n > 2, W n G 02/c(max[A/"^ fc+1 ^ 7_fc ) ra ~ <5 , 1]). Since 
7 < 1, by choosing an integer k > 1 , 1 _ 1 , we deduce that W n G Oj(l) for all n > 2 with 
/ = 2k, and we apply Lemma 4.2 to conclude. o 

4.4 Initialization 

We now establish a priori control on the correlators. We shall use: 

Lemma 4.4 Let Wn = N e for some e > 0. Assume that for any integer p, there exists 
C p > and independent of N, such that for all x G C \ [o_ , a+] .' 



/%,/3 



11/ 



dM N (Q 

x — i 



l^Nfi 



dM N (£) 



< 



C p w N 



(d(x, [o_,a+])) 



2p 



(4-9) 



T/ien, /or all n > 2, W n G 0(w%) for the norm || • || r , w/ien A^ — > 00. 
Proof. For n > 2, W n (xi, . . . , x n ) is a polynomial in: 



F;[a_,a+] 



n 






a;,-^ 



with J C {l,...,n}, and the coefficients of this polynomial are independent of N. 

Thus by Eqn. 4-9 and Holder inequality, there exists D n G 1R^_ independent of iV such 

that: 

|W n (a;i,...,a; n )| < 

( mini<i<„ d(Xi, [a-,a + \)) 

Hence, taking the sup for Xi G T, we find W n G 0(w^). 



2/i 



Lemma 4.5 Under the five assumptions of Section 4-1, Eqn. 4-9 holds for any e > 0. 

Proof. Our starting point comes from a result of Boutet de Monvel, Pastur and 
Shcherbina [dMPS95], developed by Johansson 2 [Joh98, (3.49)] and more recently in 
[KS10, (2.26)]. Let P* C Int T be a contour surronding [o_, a+\. For any : Int(r) — > C 
which is a continuous function, and real- valued on [a_,a + ], there exists a positive 
constant C such that: 



v 



1 



cxp 



2(sup z , er , \<p(z')\)w N 



^(OdJiMo - n / m&m 



<3 



2 Johansson's has written his proof in the framework [a_,0+] = M, but there is no difficulty adapting 
it to [a_ , 0+] finite. Im z should be replaced by d(z, [a_ , a+]), and its powers in the bound of his Lemma 
3. fO and 3.11 may differ, but the order of magnitude u>n (our u>n) is the same. 
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where wn = Cha.N. By Chebychev's inequality, we deduce that: 



Vt G [0,+oo[, fiYtfi 



and therefore: 



0(£)dJlMO - N / <p(OdL{0 > t( sup |0(*')l)™iv \ < 6e 



v 



VpeN, ^ 



0(OdM^(O - JV / 0(OdL(O 



<p!(sup 2 , er ,|0(z')l) 



'M^P. .,P 



W 



N 



In particular, we can apply this discussion to <j>(z) = Re -^- and (f>(z) = Im -^- where 

x is a point of I\ This leads to Eqn. 4-9. o 

In the case of a strictly convex potential, we may also use concentration of measure: 



Lemma 4.6 IfV^ is strictly convex on [a-,a + ], then Eqn. 4~9 holds with e = 0. 

Proof. Since V^°' is strictly convex on [o_,a+], V is also strictly convex on [a__,a + ] 
for N large enough. By concentration of measure, see [GZ00] or [AGZ10, Section 2.3 
and Exercise 4.4.33], there exists c > such that, for all x e C\[a_,a + ], for all e > 
and N e N: 



/%,/3 



This entails Eqn. 4-9. 



dM N (0 v 

' A*JV,/3 



dM N (£) 



> 



(d(x,[a_,a + ]))' 



< 2e~ 



4.5 Leading order of W\ 

Afterwards, all steps only rely on the analysis of loop equations. Although we already 
know the characterization of the equilibrium measure /z eq , and thus of its Stieltjes 
transform W^ , let us recall how W^~ is characterized by the loop equations. We 
write the loop equation at rank 1 (Thm. 3.1): 



AT2 



W 2 {x,x 



+{wl- 1} (x)y 



i 



(x — d-)(x — a+ 



C([a_,a+}) 2iTT X-(, (x - a_)(x - 0+) 

1 \ 1 



+ K^ lWl ) + i (l - |) (w}-%) + {x _ a _ ](x _ a+) 



N 



M vW [W\ 



{-!>! 



+ ((6- 1 W 1 )(x)) 2 -Af iSoV y[5- 1 W 1 ](x) + ^Ar iSlV y[WJ; 1} ](x 



= 
(4-10) 



2y 



We already know that the 4 th and the 5 th line are o(l). Since W 2 G o(N 2 ), W} 
satisfy the loop equation at leading order: 

^"W = (x-a-)(x-^) <4 - n) 

+ / di 1 ( ? -a-)K- a+ ) (0} ( - 1} 

Jc([a-,a+}) 2tTT X-C (x - a_)(x - 0+) 



Remark 4.1 Recall that suppyU cq = [a_,a + ] is the discontinuity locus of W{ . By 
the properties of the Stieltjes transform: 

y(x) = «) _ ^-., w 

defines a holomorphic function on U \ [a_, a+], and: 

d/ZeqfXo) 



Vx e [a_, a+], lim y(x + ze) = in 

e^o+ dx 



o 



VKe state that there exists M(x), continuous in some open neighborhood of [a_,aq_], 

y(x) = , M{x) (4-12) 

a/(x — a_)(x — a + ) 

Proof. In Eqn. 4-11, we may first deform the contour C([a_,o+]) to C([a_,a+]). 
Secondly, we can rewrite: 

(Wt 1} (x)) 2 - (vW)'(x)wl- l} (x) + ~, ^ r = 

[x — a,-){x — a + ) 

u(,) = -i+ / £ K -■»-)« -M (v/{0)yK) , v (-. lK) 

JC([a-,a+])L>C(x) zm x ~ S 

where C/(a?) is now holomorphic in some open neighborhood of [a_,a+]. So: 



y[x) 



RU) R(x) = \{x - a_)(x - a + ){V {0} )'{x) - U{x) 



(x — a_)(x — a+) ' 4 



This equation tells us that the discontinuity of y is of squareroot type. If a_ = a_ and 
a + = a_|_, we have Eqn. 4-12. If say a_ < a_, the fact that y(x) has no discontinuity 
on [a_, a_[ but a discontinuity on [a_, a + ] forces i?(x) to have a simple zero at x = a_ 
and at x = a_, so that y(x) is finite when x = a_ and vanishes as 0(y/x — a J) when 
x — >• a_. A similar statement holds if a + > a + . Then, Eqn. 4-12 holds a fortiori. o 
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4.6 First correction to W\ 

Let us reconsider Eqn. 4-10 (or the equivalent relation taking Remark 3.3 into account) 
after removing the 2 nd and the 3 rd line which has just been identified as the leading 
order. We can write as in 5 4.3: 



«+ ffi 4H)s 



<LiWi(z) = A 2 + C + D c 



where: 



8K[f]{x) = -Af {doVy [f(x)}+6- 1 W 1 (x)f(x) 



A 2 
Co 

Do 



-— 2 W 2 {x,x) 



-iri-l 



Edg T In Z 
T — 



reHard 



\Jjf (Jj — rj- t\j Ltq- 



reSoft 






By an argument similar to Eqn. 4-8, knowing that W 2 G Oi(N) implies that 
5-iWi G Oi+i(l/N). Assuming further W$ G Oy(N) implies after Section 4.3 that 
W 2 G CV +2 (1), so the 1 st line of Eqn. 4-10 is subleading compared to the 3 rd line. 
These two bounds are provided by Section 4.4 (the values of I and V do not matter 
here). Hence: 

Lemma 4.7 There exists w} 0} G H ( ^ a _ a+] such that W x = NWt ±} + w} 0} + o(l). 
Explicitly: 



w} 0} (x)=K- 1 



P 



dx 



{Wi 



{-1}, 



>)) + £ 



reHard 



a T — CL- T x — a-r 



{-lh 



+ ^(y{i>)W ')(*) 



This order was also obtained by [KS10] with similar arguments. 



4.7 Recursion hypothesis at order ko 

Let ko > — 1. We assume that the correlators W n (for all n > 1) are determined up to 
a o(iV -fe °) for the norm II ■ IL 

V / I 111 1: 



- l(k ;n) 



ko 



W n (xi,...,x n )= J2 N- k wW(x 1 ,...,x n )+N- k0 6 ko W n (xi 



, . . . , x n 



(4-13) 



k=n-2 



r{k}. 



Here, Wn (x%, . . . , x n ) are already known (they depend on f3 but not on N), and we 
call: 



to 



« 



sup ||H/W|| r 



l<ra'<n 
-Kk'<k 



l(k ;n) 
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a bound for their norm. We can always assume that l(k, n) defined for — 1 < k < k$ 
and n > 1 is an increasing function of k and n. Though the errors 8k W n are not 
supposed to be known, we assume that they are small: 

Vn>l WK^nWv <ejJ o} AJ: fe °> 

Here, e N ° depends only on N and k , and e N ° — > when N — > oo, and A„ ° is a 
constant independent of N. We may assume that A„ increases with n > 1, upon 
replacement by sup 1<n / <n A„ . When n > k + 2, we assume that Eqn. 4-13 reduces 



to: 

jxr Ar-fcn 

?fc 



w = iV- fco <j fcn w„ 



Lemma 4.4 and Section 4.7 ensure that the initial (ko = — 1) recursion hypothesis is 
satisfied. Moreover, we can take ej^ ' = 1/N, and up to a redefinition T k — > T fc _ m for 
some integer m, we can take l(— 1; n) = 4(n — 1). 

4.8 Determination of ^ W no 

Let n > 1. We now turn to the determination of the leading order of 5k W no (x,xi). 
The case (no, ko) — (1, — 1) is a bit special (because of the second term of the second 
line in Eqn. 4-10) and is given by Lemma 4.7. In all other cases, we consider the loop 
equation at rank uq (Thm. 3.2). Up to o(N~( ko ~^), the equation is true and involves 
quantities which are already known from the recursion hypothesis. The equality of 
the c^iV^ 0-1 )) involves the unknown Sk W no (x,xi). The operator /C introduced in 
§ 4.2.1 plays a special role. When the potential V has a 1/N expansion, the operator 
M introduced in § 4.2.3 also appears. We find: 

iV-to-^/C^W^foxz) = -N-** E%>\x, Xl ) - N~ k ° R^\x, Xl ) (4-14) 

with 

fco+l 

Ej*>(x, xr) := W™(x, x, x T ) - J2 V>y [w£° +1 - h} (x, x,)] 

fc=i 

JCJ fc=0 ^ 



2 y. d j Vj^fos/Ui}) Lfo)/ 1 
/3 *r~i dxi I x — X{ L(x) \x — Xi 



iei 
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and the remaining 

fc k 



i2g»>(x, xj) := 6 ko W no+1 (x, x, XJ ) + J2 N ~ k E E ^S+i^ *') ^noHJl"" >(*> x i\j) 

fc=l fe'=0 JCJ 

fco 

+ E ^'* E OW 1+0 (*» ^) ^ ji (*» ^v) 

fc=0 jcj 

+iV -fco ^ (4 o M/ |J|+1 )( x ,xj) ((JfcoW^HJiJCx.XjXj) 
jcj 

+ (l - ^^((^o^no) (»,»/) J - (8 ko W no+1 ){x,X,Xi) 

k k 

-J2 N ~ kM ( S V \ [W£H*,*l)] -J2 N ~ k ^(V^MKWn )(x,X I )} 
fc=0 ^ fc ° +i ' k=0 

2 d f ( W .,)(*,*A W ) _ XW ,_!_ + x 

p z — •* dxi x — Xi L{x) \x — Xi J 

8 2—i x - a T 

' reSoft T 

It is understood that /C and A/" g operate on the x variable. The variables xj are 
spectators. Notice that this equation is linear in Sk Q W no , up to a small quadratic term. 
Looking naively at this equation, we see that the leading term of Sk W no happens to 
be of order 1/N (giving a N~( ko+l > contribution to W no ), and is obtained by applying 
/C _1 to En ° (x,xj). To make this idea rigorous, let us bound Rn ° . Even if some 
terms in the right hand side have not been determined yet (like S). W n that we are just 
considering), we already know a bound for each of them from the recursion hypothesis. 
Very rough bounds are enough, we just need to show that the right hand side is small 
when N — > oo. When k = — 1, we must pay special attention at the terms involving 
jy-feo directly, i.e. the 3 rd line and the 7 th line in Eqn. 4-15. In the 7 th line, (5 ko+1 V)' 
is of order N^ 1 , so we obtain a term of order e N ° , which is always small. The 3 rd line 
is of order N(e N ) 2 , which is also small since we have here ejy = 1/N (Lemma 4.7). 
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For N > 1 and large enough, we have 

ll^° } llr__ < 4 fc ° } AiS 1 + iV- 1 (A:o + l)2-- 1 (* } ) 2 

,{*»} 2 no-l Ago} „£,} + (e (Jo} )2 ^-.o (A (M ); 



- i(fco + l;no) 



+£ « |l - ^ | 0< ta „> A« + AT' £ ll^ +i ,,,ll r , (> „, A™ 



fc=0 



+4 fc ° } F IMvwvllr A £° } + 4*° } N ~ k ° Wf \'llr A £° } 

A/ ^ II (.l/w; lir !(feo .„ o) n JV II ^ fco+1 yj llr i(fc ;"0) "° 



.W 2 ,... . fu, , a.. ^Pg er,^.^,!) TO I | x{/ ,,,j 

■ !(fc ;n -l) 



+% £ Q( fe0 ;no-l) ^|c| + Cto-D+i j^^^|X^)l I A,- 1 

I 2 7n#50ft fc0+n , 1 ^^ 

/3 d{T,[a-,b-]) n 

Given the control provided by the recursion hypothesis, this inequality is correct pro- 
vided we choose: 

l(k + 1; n ) > max [l(k ; n - 1) + 2, /(£; ; n ) + 1, Z(fc ; ^o + 1)] (4-15) 

Accordingly, R { ^ } ->■ when JV -> oo. Eqn. 4-14 tells us that E^ o} + Rn o} G Im/C 
for any N. Since Im/C is closed (Lemma 4.1), we know that En ° G Im/C, and also by 
difference Rn ° G Im/C for any N. And, by continuity of /C _1 , we deduce: 

h W no = ^ n { fco+1} + 1 4o + i^„ 

where: 

^o + i} = -/C-^oi], 4 o+1 ^ no = -XT 1 ^] G o(l) (4-16) 



The previous inequality is more precise about the o/(fc 0+ i ;no )(l): there exists a constant 
A£ 0+1 \ such that 

ll^ + i^ || rufco+i;no) <ASo + i} max(iV -i ;e ifo }) 

4.9 Remarks 

The recursion hypothesis tells us that W^o = whenever n > ko + 2 (we call *[ko) 
this recursive assumption). Let us see what happens at order fc + 1 (here, ko is fixed, 
but no is free), by looking at Eqn. 4-16. 

• The term W^+i vanishes whenever no > k + I. 

• The term Wj,^ W^ \J may be non zero in case k + 1 > \J\ > uq — ko — 2 + k. 
This is impossible to fulfil as soon as n > k Q + 3. 
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The term ( Wn ) vanishes whenever no > ko + 2. 

The term involving W„ _■ 1 vanishes whenever no > k$ + 3. 



n - 

Accordingly, w£ 0+1} = when n > (k + 1) + 2, i.e. *[k + 1] holds. This is 
just the manifestation of Lemma 4.2. Hence, we have propagated the full recursion 
hypothesis to order k^ + 1 . An easy recursion shows that Wn are actually holomorphic 
functions on the domain C\ [a_, a+], i.e belongs to the subspace 'H n ., a _ a i of % n .\ a _ a y 
Therefore, we can contract the contour to C([a_, a+]) in the expression of /C _1 (Eqn. 4- 
5) when computing Wn with formula 4-16. 

Since l(k;n) = 4(n — 1), the minimal solution of Eqn. 4-15 is l(k\n) = 4(n + k). 
Indeed, in this proof, we need to have a more restrictive control on the error done at 
height n + k, in order to bound the error done at height n + k + 1. Nevertheless, since 
T^ C Int(r^) for all /, we can at the end make the weaker statement that, for any n 
and k: 

\\S k W n \\ TE -+ (4-17) 

when A^ — > oo. However, we necessarily have d(Ti,Ti +1 ) — > when / — >■ oo, so that 
the constant Q which allows us to bound the derivative of a function with the function 
itself (Eqn. 4-3), blows up. This means that Eqn. 4-17 cannot be uniform 3 in n and k, 
even when (3 = 2. 

A posteriori, from Eqn. 4-17, we can deduce by choosing rather l(k ; n ) = 8(n +k ): 



\5. 



Subsequently, upon redefinition of the constant A„ ° , we may choose e N ° = l/N. 
Finally, we can make the weaker statement that, for any n and k: 

\\S k W n \\ TE e o(l/N) 

without uniformity in n and k. 

5 Proof of the main results 

5.1 Expansion of the correlators 

We wish to study the /3 ensembles on a given interval [&_, b+], with the hypotheses 1.1 
on the potential V. When both edges are hard, Hyp. 1.1 are equivalent to the five 
assumptions of Section 4, so the Proposition 1.1 is already proved, as we have shown 



3 We thank Pavel Bleher for pointing out a mistake in a former version of the article, which we 
corrected by introducing this family of nested contours. 
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recursively that Eqn. 4-13 holds for all ko- Let us now assume that one of the edge is 
soft. The equilibrium measure /i eq := /i eq "' with support [a_,a + ] C [&_,&+] also 
coincides with /x eq , where a_ can be any point in [&_, cn_[ if b_ is a soft edge, and 

a_ = 6_ else (resp. a + can be any point in ]a + , b + ] if b + is a soft edge, and a + = b + else). 
When b T is a soft edge, "offcriticality" implies that S(x) is positive in a neighborhood 
of a T in [6_, 6 + ]\]a_, a+[. So, one can choose an interval [o_,04-] C U, and such that 
the five assumptions of Section 4 are satisfied for dfi N o~' a+ . Theorem 4.1 then can be 
applied: there exists an asymptotic expansion 

Wr' [a -' a+] (x u ...,x n )= J2 N- k W^ {k} (x l ,...,x n ) (5-1) 

k>n-2 

with respect to the norm || • || r where Te Q U can be any contour surrounding 
[a_,a + ] but not the zeroes of S. The "large deviation control" on [&_,&+] allows to 
use Proposition 2.2: there exists rj > such that, for any contour T' E C C surrounding 
[&_, b + ], there exists T n<r > such that: 



| W ,V;[6_, 6+ ]_ W ,V;[a_,a + ]|| < j , fi - 



iV?j 



This implies that the right hand side of Eqn. 5-1 is an asymptotic series for 
W n ~' (a?i, . . . , cc n ), uniformly for (a?i, . . . , x n ) in any compact of (C \ [6_, 6+])". 

We give below a more transparent condition, which imply the "large deviation 
control" assumption on [&_,&+]: 

Remark 5.1 If S(x) > whenever x E [6_,6 + ], i/ien j7" y; ' 6_ ' 6+ ^ achieves its minimum 
value only on [a_,a+], 

Indeed, l 7" V; [ 6_,6+ 'l(x) is differentiable when x G]6_,6 + [\[a_,a+], and we have: 

(^[ b -' b+ ](^))' = {Vi ° } J {x) - Wt 1} (x) = y(x) = S(x)a(x) 

The sign of the square root a(x) is determined for example by the positivity conditions 
1-2 on J~ v n b -' b +i. If we assume that S do not vanish on [6_,6 + ], this implies that 
jv-,[b-,b + ] j g s trictly decreasing in [6_, a_[ and strictly increasing on ]a + , b + ], hence the 
remark. 

5.2 Expansion of the free energy 

So far, we only have determined the expansion of the correlators which are by definition 
derivatives of the free energy. To find the free energy itself, one would like to interpolate 
between our initial potential V, and a simpler situation, using that the difference 
depends on the correlators. For any fixed a_ < a + , and fixed nature of the edges 
X-t G {hard, soft}, we denote by V^' x ^ the set of potentials V: 
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• defined at least on some interval [o_,a+] 5 [a_,a+], with a T ^ a T if X T = soft, 
and a T = a T if X T = hard ; 

• which satisfies the five assumptions of Section 4.1 on [a_,a+], in particular is 
offcritical on [a_,o+] ; 

• for which the equilibrium measure fJ, e q has [a_,aq_] as support, 

• and such that a T is an edge of nature X T . 

Lemma 5.1 V" + ' Y + is a convex set. 

Proof. Let V , V x G V°+{+, and set V, = (l-s)Vb + sVi for s e [0, 1]. V and V 1 are at 
least defined on a common interval [o_, a+] 2 [a_, a+]. Let us call v s = d/i e q the 

equilibrium measure for the potential V s on [o_, o + ]. We observe that (1 — s)dz/ + sdv\ 
is a probability measure which is solution of the characterization of dL s by Thm. 1.1. 
Therefore, dL s = (1 — s)dL + sdL±. Besides, we know that there exists a function S s , 
regular in a neighborhood of [a_,a + ] in the complex plane, positive on [o_,a+], such 
that: 

for s = or s = 1. Since the edges are of the same nature in Vo et VI, we must 
have S s — (1 — s)So + sSi. Since So and Si are positive on [a_,a_|_], so is S s . Hence 

Corollary 5.1 Lei Vo, V\ G V"^' x *. When a_ ared a + satisfy the condition above, the 
quantity: 

lnZ V^-,a +] _ lnZ Vof-,a + ] = _^l f\l ^_ (Vi(0 _ Fo(e))w K ; [a-, a+ ] (e) 

2 JO •/c([o_,o+]) 2* 71 " 

/tas a /arge iV asymptotic expansion of the form: 

\nZ^ - InZ^p = J2 N~ k F^ Vv ' [a -' a+]m 

k>-2 

where: 

1 fc+2 

F jo-^vi ! [«-,a + ] ! {fc} = _/? I" d / d^ y> (yW(e)-V^ m} (0) (^^ a - a +]) {fc+1 - m} (^) 

2 Jo /C([ a _,a + ])2Z7T^ 
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Proof. Since V s satisfies the five assumptions of Section 4.1 for any s € [0,1], we 
can apply our main theorem to W 1 s . Moreover, since we do not reach a critical point 
when s is in the compact [0, 1], we know that the error 0(N~ K ) made if we replace 
W 1 " by X^fcJli N~ k W 1 s is uniformly bounded with respect to s on some contour 
surrounding [a_,04_] and in the analyticity domain of V. Therefore, we can exchange 
the integral and the sum in the asymptotic expansion. o 

For instance, when V satisfies the five assumptions of Section 4.1 on some interval 
[a_,a + ], such that a± are soft edges, one can interpolate between V and a gaussian 
potential corresponding to an equilibrium measure with support [a_,a+]: 

I r — 



(a + — a_ 



Proposition 5.1 Let V be a potential satisfying the five assumptions of Section 4-1 on 
some interval [a_,a + ], such that a± are soft edges. For all s G [0,1], (1 — s)V+sVQ ta _ ta+ 
belongs to V"+s° ft and we have the following asymptotic expansion when N — > oo; 



^,, = <r ' a+ -p ( £ at-f; 



,V __ ry V G,a_,a + ^^ ( ST^ A] -^k T7 V ^ V G,a_ , a+ ;[a~ ,a + Y,{k} 

' fe>-2 



where the pref actor is a partition function of the gaussian (3 ensemble (see Eqn. 1-4): 



Z 7V,/3 - Z iV,G/3E ^ j 

According to the discussion of § 5.1, we can weaken the hypothesis of the proposition 
above to find Theorem 1.2. 

5.3 Central limit theorem 

Eventually, our results imply the central limit theorem proved by Johansson [Joh98], 
but here integration is taken on a compact set [o_, a+] instead of the real line (in fact 
as our derivation is quite similar to Johansson's, this is not surprising). For simplicity, 
we take here the hypotheses of Section 4, although we could refine to hypotheses 1.1 
following § 5.1. 

Let h : [o_, o+] -4 1 be a function which can be extended as a holomorphic 
function defined on some neighborhood of [a_, o+], let us take V = V^ independent 
of JV, yW = | h and define V h = V™ + iV" 1 ^ {1} = V-^h. Then: 



y ; [a_,a+] 



N 

exp ( £ h(Xi 

i=l 



yV h \[a-,a+\ 
7 V;[a_,a+] 
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and we can use Corollary 5.1 to derive its large N asymptotics. Indeed, we have 

N 



ln/x^ 



V;[a-,a+] 



exp 



I>(A< 



d.s 







^ rtrV. 



C({a-,a+}) 



2m 



w^iOHO 



By Theorem 5.1, or simply at the point of Lemma 4.7, we have: 

d/icq(^) 






€-v 



dx 



(wl 



V;{-1], 



X ) + > 



reHard 



\Jbrr- tX -j- Ju (X<r- 



2,s 



- f^w?*- 1 ^*) 



< 



/vw; 



V sh ;{-1} 



W v eh -M + o(1) 



which shows the: 

Proposition 5.2 Central limit theorem. 



! y;fa_,aj 



N 



cxp 



J>(Ai) 



i=l 



1 



N / d/ie^r?) %) + m[h] + - C[/i] + o(l) 



to£/i m[/i] t/ie linear in the function h, given by: 



m[h] 



1 



/3 



C([o_,a+]) 



2lTX 



K -.{A K m-., w)+e _^ 



reHard 



(X-j- (X -T- tVy (X-t 



HO 



and C[h] the quadratic function of h given by 

2 



C[h] 



P 



C([a-,a+}) 



2«7T 



^;{-l} 



A4W ;1 ') (0M0 



Therefore Y2 i=1 h(\i) — N f d/j, eq (r)) /i(r/) converges towards a Gaussian variable with 
mean m[h] and covariance C[h]. 
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